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THE MINIMAL ANGLE CONDITION FOR QUADRILATERAL FINITE 
ELEMENTS OF ARBITRARY DEGREE 

GABRIEL ACOSTA AND GABRIEL MONZON 


Abstract. We study Lagrange interpolation error estimates for general quadrilateral Qk 
finite elements with k > 2. For the most standard case of p = 2 it turns out that the constant 
C involved in the error estimate can be bounded in terms of the minimal interior angle of 
the quadrilateral. Moreover, the same holds for any p in the range 1 < p < 3. On the other 
hand, for 3 < p we show that C also depends on the maximal interior angle. We provide some 
counterexamples showing that our results are sharp. 


1. Introduction 

This paper deals with error estimates in the norm for the Qk Lagrange interpolation on 
a general convex quadrilateral K C IR^. Denoting the interpolant with Qk the standard error 
estimate is nsnally fonnd in the form 

||tt Qfc^llo,p,A h\u p p' ^ Ch Ittl(I'l) 

being h the diameter of K. Ineqnality (|l.ll) involves the error estimate 

||tt Qk^tWop}^ ^ Ch Iu] 

and the seminorm estimate 

Qkulip j^ ^ Ch (1'3) 

A central matter of (HI]) concerns the dependence of C on basic geometric quantities of the 
nnderlying element K. It is known that the constant C in (|1.2I) remains nniformly bounded 
for arbitrary convex qnadrilaterals (see Theorem EH). However this statement is false for the 
constant C in (ll.^p (see for instance the connterexamples in the last section). The primary goal 
of this paper is to stndy the dependance of C in (|1.3p on the interior angles of K. Althongh the 
role of the interior angles have been related to C in many previons works, none of them, to the 
best of the anthors knowledge, have given a result as plain as the one offered in this paper. For 
instance, bonnding the minimal and the maximal inner angle is considered a central matter in 
mesh generation algorithms since the early work by Ciarlet and Raviart [8], however no proof 
of snfficiency has been given as far (at least for an arbitrary degree of interpolation). 

In order to present onr results let ns first introduce the following classical definition that we 
write for both triangles and qnadrilaterals for further convenience. 

Definition 1.1. Let K (resp. T ) be a convex quadrilateral (resp. a triangle). We say that K 
(resp. T) satisfies the minimum angle condition with constant G IR? or shortly mac{fi)rn), if 
for any internal angle 9 of K (resp. T) ^ 9. 
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Our first result says that the constant in (11.31) . for a fixed degree k and a fixed value of p 
with 1 < p < 3, can be written as C = C{'ijjm)- As a consequence, the same can be stated 
about the constant in (jl.ip . This seems to be the most general result available for quadrilaterals 
in the case k > 2. In spite of the fact that much weaker geometrical conditions are known to 
be sufficient for the case A: = 1, we show, by means of a counterexample, that they fail for a 
higher degree interpolation. This counterexample also warns that removing the minimal angle 
condition may indeed lead to a blowing-up constant in (II.3p . 

The mac is the most standard condition considered in textbooks for triangular finite elements. 
Actually, in that case, it is equivalent to the so called regularity condition, i.e. equivalent to the 
existence of a constant a such that 

h/p < a, (1.4) 

where p denotes the diameter of the maximum circle contained in T. On the other hand, the 
term anisotropic or narrow is usually applied to elements that do not satisfy (ILI. Even when 
triangles can become narrow only if the minimal angle is approaching zero a very different 
situation occurs on quadrilaterals. Indeed, in that case the mac condition is less restrictive than 
(jl.4l) since arbitrarily narrow elements are allowed with a positive uniform bound for the minimal 
angle (for example, anisotropic rectangles always verify mac{'K/2)). Anisotropic elements are 
important for instance in problems involving singular layers and the first works dealing with 
them arise during the seventies showing that (11.41) can be replaced (for triangles) by the weaker 
following condition 

Definition 1.2. Let K (resp. T) be a convex quadrilateral (resp. triangle), we say that K (resp. 
T) satisfies the maximum angle condition with constant G IFl, or shortly MAC{fiM), if for 
any internal angle 6 of K (resp. T) 6 < < t. 

Indeed, in [3 [9] it is proved that the MAC is sufficient to have optimal order error estimates 
for Lagrange interpolation on triangles. In the case of quadrilateral elements, m it is also a 
sufficient condition as it was shown by Jamet m for k = 1 and p = 2. This condition is less 
restrictive than that propoposed in [8] where the authors require the existence of two positive 
constants pi,p 2 such that 

h/s < Pi (1.5) 

where s is the length of the shortest side of K, and 

I cos(0)| < P 2 <1 (1-6) 

for each inner angle 9 of K. Observe that under the regularity condition (|1.4p the quadrilateral 
can degenerate into a triangle (for instance if the shortest side tends to zero faster than their 
neighboring sides or if the maximum angle of the element approaches vr), however this kind of 
quadrilateral cannot become too narrow. Condition (II.6p will play an important role in the 
sequel and therefore we introduce the following alternative definition. 

Definition 1.3. We say that a quadrilateral K satisfies the double angle condition with con¬ 
stants if rn,'<pM, or shortly DACif K verifies mac{'iljm) andMAC{'fiM) simultaneously, 
i.e., if all inner angles 6 of K verify 0 < < 9 < < "t. 

The DAC allows anisotropic quadrilaterals (such as narrow rectangles) as well as families of 
quadrilaterals that may degenerate into triangles. To see that consider, for instance, a quadri¬ 
lateral with vertices (0,0), (1,0), (s, 1 — s) and (0,1 — s) and take 0 < s —>■ 0. 

For anisotropic quadrilaterals several papers have been written mainly in the isoparametric 
case with A; = 1. In [nmi] narrow quadrilaterals are studied and the authors require the two 
longest sides of the element to be opposite and almost parallel, the constant C obtained by them 
depends on an angle which in some cases is the minimum angle of the element. Anisotropic 
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error estimates for small perturbations of rectangles have been derived in EE]. On the other 
hand, for k = 1, more general and subtle conditions can be found in the literature. For k = 1 
and p = 2, it is proved in [2] that the optimal error estimate (11.311 can be obtained under the 
following weak condition 

Definition 1.4. Let K be a convex quadrilateral, and let di and d 2 be the diagonals of K. We 
say that K satisfies the regular decomposition property with constants £ IR and 0 < fjM < t, 
or shortly RDP{N,’iPm), if we can divide K into two triangles along one of its diagonals, that 
will be called always di, in such a way that |(i 2 |/|di| < N and both triangles have its maximum 
angle bounded by ^m- 

In Remarks 2.3 - 2.7 of [2] it is shown that the regular decomposition property RDP is 
certainly much weaker than those considered in previous works (including [uiiiEKinKiaiii]). 
We collect for further reference some elementary remarks 

Remark 1.1. If a quadrilateral K satisfies (11.41) then K verifies RDP{a,'ijj) whith = fiicr). 
Indeed, m implies that K verifies maciO) for some 0 = 0[a) > 0. Therefore there is at most 
one angle of K not bounded by ir — 9. Dividing K by the diagonal di containing the vertex 
associated to that angle we get that K satisfies RDP{a,7r — 9). 

Remark 1.2. If a quadrilateral K satisfies MAC{'ijjM) then K verifies RDP{l,fjM), cis one 
can see by taking di as the longest diagonal of K. 

Since, by definition, DAC{firn,'4’M) implies MAC{'fiM) we have 

Remark 1.3. If a quadrilateral K satisfies DAC{'iprn-,'f'M) then K verifies RDP{1 ,'iIjm)- 

Remark 1.4. If K verifies the mac{'il)m), then K either verifies 

(1) DAC{iim,T^ - ^) or 

(2) the regularity condition (II. 4p with C = C{firn)- 

Indeed, assume that (1) does not hold. Then K has an internal angle which is greater than 
T — it is easy to see that this angle is unique so we can call it 9. Divide K into two 

triangles Ti and T 2 through the diagonal opposite to 9 in such a way that 9 becomes an internal 
angle ofTi. Calling fii and /I 2 to the other angles ofTi it follows that fii < 'firn/‘2‘ with i = 1,2. 
Let i = 1,2, the complementary angle of fii (w.r.t. the corresponding internal angle of K). 
It is easy to see that 7 * > 'fiml‘2‘, so T 2 is a triangle that have its three internal angles bounded 
away from 0. To be more precise, T 2 verifies mac{'il)rn/‘2), therefore T 2 is a regular triangle in the 
sense of (HI with C = C'(V’m/2). From this fact follows easily that K is a regular quadrilateral 
in the same sense, i.e., in such a way that o holds (actually the same constant C can be 
used). □ 

Remark 1.5. Combining Remarks II.il 11.31 and \1.4\ it is clear that mac RDP. 

An strikingly result is that the RDP, in spite of being appropriate for p = 2, does not work 
for arbitrary values of p. Indeed in [3] and for k = 1 the results of [2] are extended for the error 
in with 1 < p < 3. Moreover, it is shown, by means of a counterexample, that this range 
for p is sharp. As a consequence the stronger DAC (i.e. (ll.Op 'l is proposed and shown that 
under this condition the error estimate hold for all p > 3. The second resnlt given in the present 
paper is that, for p in this range, the DAC is also a sufficient for any k >2. 

For the reader’s convenience we summarize in the following table simple and sufficient geo¬ 
metric conditions pointing out the role of p and k 



I < p < 3 

3 < p 

k = 1 

RDP 

DAC 

k>2 

mac 

DAC 







THE MINIMAL ANGLE CONDITION FOR QUADRILATERAL FINITE ELEMENTS 


4 


the first row of the table was proved in [3] while the new results are given in the second row. 

To finish this short review we recall that for k = 1 more results are available. In 
error estimates are obtained for the Qi isoparametric Lagrange interpolation under a weaker 
condition than the RDP. This condition can be regarded as a generalization of the last one and 
therefore called GRDP. 

This paper is structured as follows: in Section 2 we introduce a family of reference elements 
appropriate for dealing with general convex quadrilaterals and some key results are provided. 
In Section 3 and Section 4 our family is related to different geometric conditions (such as RDP, 
MAC, DAC and mac) while some properties about the distribution of the interpolation nodes 
of the family are studied. Section 5 gives the general approach for bounding the interpolation 
error. Finally the main results as well as the counterexamples can be found in the last section 
of the paper. 


2. The family of reference K{a,b,a,h) 

With K we denote an arbitrary convex quadrilateral with vertices Vi,V 2 , 1 / 3,14 enumerated 
in counterclockwise order. For positive numbers a,b,a,b, we use K{a,b,a,b) to represent a 
quadrilateral (always convex) with vertices I 4 = (0,0), 14 = («, 0), V 3 = (a, 6) and 14 = (0,6). 
In particular, K = K{1,1,1,1) is the reference unit square and for any positive integer k we 
consider {k + 1)^ points of coordinates x = j/k and y = i/k. For K we write 

Vi = Moo, 14 = Mofc, 14 = Mkk and I 4 = Mko- ^ ^ 

We define as usual Fk : K ^ K as Fk{x) = Yli=i ^i4>iix) being (pi the bilinear basis function 
associated with the vertex Vi, i.e., <pi{Vj) = <5^. Then, in K we have {k-\- Vf points {Mjj}o<ij<fc 
defined by 

Alij = FxiMij). 

For quadrilateral elements (isoparametric when k = 1 and subparametric otherwise), we 
have the basis function on K defined by (l)ij{X) = (pij{F^^(X)) where <pij € Qk{K) verifies 
(f)ij{Mir) = and therefore the Qk interpolation operator Qk, on K is given by 

Qku{X) = QkHX) 

where X = Fk{X) and Qk is the Lagrange interpolation of order k oi u = u o Fp^ on K. 
Interpolation nodes of the form {Mjj}i<jj<fc_i are called interior and any Mij which is not 
an interior node is called an edge node. Also of interest is the triangle T{a,b) of vertices 
= (0,0), 14 ^ = (a, 0), = (0, 6). The interpolation nodes of the Lagrange interpolation 

operator 11^ € Pfc of degree k are given by = [aj jk, bi/k), D < i + j < k. With C we denote 
a positive constant that may change from line to line. Sometimes we also use the notation x ~ y 

for positive variables if they are comparable in the following sense < y < Cx. 

For any element of the type K[a,b,d,b) considered in this work the following condition will 
become relevant in different contexts 

Condition (AI): —,t<C. (2.1) 

a b 

This condition takes sometimes the more restrictive form 

Condition (Dl): -,7<1. (2.2) 

a b 

In spite of the fact that both (12.11) and (12.2p look similar they characterize, under the supple¬ 
mentary geometric assumption (12.3p (see below), different classes of elements. 
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Calling di to the diagonal joining V 2 and V 4 we see that di divides K into two triangles, that 
we call Ti and T 2 (see Figured]). For the angle a of Ti placed at V 4 , we introduce 

Condition {D2)-. —< C, (2.3) 

sin{a) 

which says that a is bounded away from 0 and tt. 



Figure 1. Notation for an element K{a,b,d,b). 

Finally, in order to exploit some results given in previous works, we introduce yet another 
useful condition 

Condition (A 2 ): |/| < Clsl. (2-4) 

where I is the segment V 3 V 4 joining V 3 and V 4 and s denotes the shortest side of K{a,b,d,b). 
That is, (A2) amounts to say that the side I is comparable to the shortest side of K. 

Not difficult to prove is the following 

Lemma 2.1. Let K{a,b,d,b) a general convex quadrilateral. Then, conditions [A1,Z12] and 
[A2, D2] are equivalents. 

Proof. That (A2) implies (Al) follows easily as we have a < |/| < Clsl < Ca and similarly 
\b — b\ < Cb. Hence triangular inequality yields b < Cb. On the other hand assume [Al,H2]. 
Thanks to {D2), and using the law of sines for the triangle A(V 2 , V 3 , V 4 ) we see that the angle 
/3 at V 3 is away from zero and tt. Indeed < C due to (Al). Since both a and 

/3 are away from zero and tt the law of sines says that actually di is comparable to IV 2 V 3 I. It 
implies, in turn, that |/| < Cmin{|V 2 F 3 I, Now consider two cases: if ^ approaches zero 

then necessarily > C > 0, otherwise the angle a can not obey (D2). Hence we have 

simultaneously, b < a and (due to (Al)) |/| = \J {h — hY + a? < C\b — b\ < Cb showing (A2). To 

finish, let us assume that - > C > 0 , in such a case if a and b are comparable we have nothing 
to prove since then they are also comparable to di. Therefore we may assume that | approaches 

cero. In that case < C since a can not get close to vr nor to zero. Therefore we have a < b 

and |/| = \l{h — b^ + a? < Cd < Ca thanks to {DI). The lemma follows. □ 

The class of reference elements of the type K{a,b,d,b) is adequate for dealing with several 
geometrical conditions. For instance, as we show below, [Al,Z12] (or equivalently [A2,D2]) 
describe any element under the RDP and [HI, D2] the family of elements under the DAC, etc. 
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In order to do so, the key tool is given by the following elementary lemma 

Lemma 2.2. Let K, K he two convex quadrilateral elements, and let L : K ^ K be an affine 
transformation L{X) = BX + P. Assume that L{K) = K, ||i?||, < C (in particular the 

condition number k{B) < C). If Qf. is the Qk (isoparametric for k = 1 or subparametric for 
k > 1) interpolation on K and u = uo L~^ then for any p > 1 

Ci\u- QkU\l^pJ< < \u- Qku\l,p,K < ^ 2111 - QkU\l,p,K 

and 

— W\m,p,K < ^*2 

for any m > 1 . 

Proof. By definition of the interpolation we have 

Qku{X) = Qku{Ff,\X)) and Qja 1{X) = Qkil{F^\X)). 

Where X denotes the variable on K. Since L is an affine transformation, F-j^ = L o Fk and 
so QiM{X) = Qku{X). Then, the lemma follows easily by observing that ||i?|| , ||i?“^|| < C. □ 

Definition 2.1. We say that two quadrilateral elements K,K are C-equivalents (or simply 
equivalents ) if and only if they can be mapped to each other by means of an affine transformation 
of the kind given in Lemma \2.2[ 

Taking into account that each geometric condition defined as far is going to be mapped to 
an appropriate class of equivalent elements K{a,b,a,b) it is important to consider the map 
Fk '■ K —>• K{a,b,a,b), 

Fk{x, y) = (ax(l - y) + axy, by{l - x) + bxy) = {x, y), (2.5) 

as well as its associated Jacobian 


DFK{x,y) 


a + y{a — a) x(d — a) \ 
y(b — b) b + x{b — b) ) ’ 


( 2 . 6 ) 


Jk := det(DFK)(x, y) = ab{l + x{b/b - 1 ) + y{d/a - 1 )). (2.7) 

Observe that since K is convex, we have Jk > 0 in the interior of K. Indeed, since Jk is an 
affine function it is enough to verify that it is positive at some vertex of K and non negative 
at the remaining ones. The positivity at Vi = (0, 0) is trivial, as well as the non negativity at 
1^2 = («) 0) and Vi = (0, b). On the other hand, since K is convex, (a, h) lies above the segment 
joining V 2 and V4 (for this segment y(x) = —{b/a){x — a) and 0 < ^ ~ I + f “ therefore, 

Jx(l, 1) = ab{d/a + b/b — 1) > 0. (2-8) 


Following again mm we introduce for p > 1 the next expression that becomes useful in the 
sequel 

Ip = Ip{a,b,d,b) := [ [ - — ^ -- —- dxdy, (2.9) 


0 Jo {1 + x{b/b - 1) + y{d/a - 1))P ^ 
where the numbers a, b, d, b are compatible with an element K{a, b, d, b). 


Lemma 2.3. If K = K{a,b,d,b) is convex and (Al) given by (12.11) holds, then for any p > 1 
and for any function basis f there exists a positive constant C such that 

df ^ „ b 

0,p,K 


dx 


< C- 


aP 


-1 


( 2 . 10 ) 
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d(f) 


dy 


< C 


h. 


0,p,K 


Jjp-1 P 


( 2 . 11 ) 


Proof. Let (p the function basis on K corresponding to p, then (from the chain rule) follows 
that 




1 


JK{x,y) 


b + x(b — b) —y(b — b) 
—x{a — a) a + y{a — a) 


( ^ 


where Jk{x, y) = det(IlF/^) = ab{l + x{b/b — 1) + y{a/a — 1)). 


Calling R{x,y) = (1 + x{b/b - l))#(x,y) - y{b/b - l)||(x,y) and S{x,y) = -x{a/a - 


d^. 


dy' 


l))§|(®,y) + (l + y(a/a- l))|f(^,y) we have 


dp I 


dp 

dx 


O Fk (x, y) = 


JK{x,y) 


R{x, y) and 


o Fk {x, y) = 


By changing variables we get 
dp b 


and 


dx 

dp 


1 /■! 


dy 


0,p,K 

P 

0,p,K 


dp 
dy 


V/o Jo [I + x{b/b - 1) + y{a/a - l))P 

\S{x,yr 


a 


JKix,y) 

— dxdy 
-1 " 

dxdy, 


S{x,y). 


Jo Jo {I + x{b/b - 1) + y{d/a - 1))P-^ 
and the proof concludes using that R and S are uniformly bounded since they are polynomials, 
0 < X, y < 1 and 0 < d/a, b/b < C by (Al). □ 

Previous result provides bounds for any basis function. As we show later basis functions 
associated to internal nodes require a particular treatment. In particular we have, 

Lemma 2.4. If K = K{a,b,d,b) is convex and (Al) given by (12.Ih holds then for any p >1 
and for any function basis p associated to an internal node of K, there exists a positive constant 
C such that 

dfh P hr- _ 1 

<C—r \1 — b/b\^Ip +max{l,{b/b)^~^/2} 


dx 


dp ^ 


dy 


< C 


aP-i 
a 


0,p,K 


bp 


-1 


(d/ayip + max{l, {b/by ^/2} 


( 2 . 12 ) 

(2.13) 


Proof. Since p is associated to an internal node on K, it follows that p is associated to an 
internal node on K so that there exists P G Qk- 2 {K) such that p{x, y) = x(l —x)y(l —y)P(x, y). 
Then 


dp 

dx 


= y{l — y)A and 


dp 


= x(l — x)B 


where A{x, y) = -^ [x(l - x)P(x, y)] and B{x, y) = ^ [y(l - y)P{x, 


From the chain rule follows that 

1 


( (dp jp \ 

gpoFK (x,y) 




JK{x,y) 


dy 

b + x(b — b) —y(b — b) 

—x(a — a) a + y{d — a) 


y(i - y)Mx,y) 

x(l — x)B{x, y) 
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where Jk{x, y) = det(-DF/^) = ab{l + x{b/b — 1) + y{a/a — 1 )). 
Calling S = xy[{l — x)B — (1 — y)A] and R = yA we have 

d(f) 
dx 


o Fk ) (x, y) = (1 - + (1 - y)R 


and by a change of variables we get 


d(f 


{l-b/b)S + {l-y)R 

p 

dx 

0,p,K Jo Jo oP ^ (1 + 

x(b/b - 1 ) + y{d/a - l))P“i 


Using the fact that S and R are uniformly bounded we see that 


d4> 


dx 


1 /■! 


< C 


0,p,K 



\l-b/b\P + {l-y)P 


Then 


d(j) 


dx 


< C- 


0,p,K 


oP 


-1 


./o Jo (1 + x{b/b - 1 ) + y{a/a - 1 ))p-i 

{i-vY 


dxdy. 


dxdy. 


1 


\l-b/b\Plp + 



— dxdy 
0-1 ^ 


/o Jo {1 + x{b/b - 1) + y{d/a - 1))P 
From the convexity of K we have a/a + 6/5 — 1>0, hence 

1 + x{b/b — 1) + y{d/a — 1) > 1 + x{b/b — 1 ) — yb/b 
Assume now that b/b < 1. Since 0 < x < 1 we conclude 1 + x{b/b — 1) > b/b and finally 

1 + x(b/b - 1 ) + y{d/a - 1 ) > b/b{l - y). 

Therefore 

-1 2 V 6 


1 /■! 



'o Jo {l + x{b/b-l) + y{d/a-l))P 
On the other hand, if 6/6 > 1 


hence 


1 + x(b/b — 1 ) + y{d/a — 1 ) > 1 + yidja — 1) >1 — y 

(1 - yY 


1 /■! 


0 Jo {l + x{b/b-l) + y{d/a-l)Y 


-1 


dxdy < 1 , 


and (I2.12P follows. 

Finally, the estimate for 


90 

dy 


0,p,K 


can be obtained in a similar way from the expression 


(if ^ ^ -y)^] 

where R = x(l — x)B + xR. □ 

In order to clarify in advance the role of the term | in the previous lemma let us notice the 
following 


Lemma 2.5. For any arbitrary convex quadrilateral K(a, 6 , a, 6 ) under condition [Al, D2] (equiv. 
[A2, D2]) there exists another equivalent element (in the sense of Definition \2. 1\) obeying [Al, D2] 
(equiv. [A2,L)2], see Lemma \2B\) with the same constants and for which 5 > 
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Proof. Consider the triangle V 2 V 3 V 4 and the angles a and /3 at V 4 and V 2 respectively. If for the 
original K{a,b,d,b), ^ then we have nothing to prove. Otherwise \ < \ and hence we see 
that a < (5. On the other hand, since both are interior angles of a triangle, /3 < vr — a, therefore 
using {D2) we see that (3 is away from 0 and vr and therefore under a rigid movement we can 
transform our element into K(b,a,b,d). The resulting element satisfies the required conditions 
[A1,T)2] with the same constants than those K(a,b,d,b) and the lemma follows thanks to the 
fact that f > 5 (see () 2 . 8 p l. □ 

Remark 2.1. Observe that previous lemma also applies to elements K{a, b, d, b) under [Dl, D2] 
since Dl Al. 


3. K{a,b,d,b) AND different geometric conditions 

In this section we explore in detail how to use the family K{a, b, d, b). The following lemma 
is useful in the rest of this section. 

Lemma 3.1. Let L be the linear transformation associated with a matrix B. Given two vectors 
wi and W 2 , let a be the angle between them and let a be the angle between L{wi) and L{w 2 )- 
Calling k{B) the condition number of B we have 

2 f 2 \ 2 

k[B)tt \ K[B)7r J k[B)tt 

Proof. The proof is elementary and can be found in [T]. □ 

3.1. The RDP and the family K{a,b,d,b). In order to characterize the elements under the 
RDP we begin with the following elementary result 

Lemma 3.2. Let K be an element of the type K{a,b,d,b) and assume [A 1,112] (equivalently 
[A2, D2]). Then K verifies the RDP with constants depending only on those given in conditions 
[A1,D2]. 

Proof. Follows straightforwardly taking di = V 2 V 4 as the dividing diagonal. □ 

Now we are ready for the following characterization 

Theorem 3.1. Let K be a general convex quadrilateral. Then K verifies the RDP if and only 
if K is equivalent to some K{a,b,d,b) under [A2,il2] (equiv. [A1,D2]). 

Proof. First we assume that K is equivalent to some K{a, b, d, b) under [A2, D2]. From Lemma 
we know that K{a,b,d,b) verifies RDP{N,'ijjM) with constants bounded in terms of those 
given in [A2,Z12]. Since K = L{K{a,b,d,b)) for certain affine mapping Lx = Bx + P of the 
type considered in the Dehnition 12 . 1 1 we see from Lemma l3. II and taking into account that such 
an L preserves lengths (up to a constant depending on ||P||, ||P~^|| < C) that K verihes the 
RDP with constants depending on L as well as the RDP constants associated to K{a,b,d,b). 
To show the other implication we follow [2]. Assume that K satisfies the RDP and divide it 
along dl into Ti and T 2 in such a way that all their interior angles are bounded by while 
< N. We choose the notation in such a way that the shortest side of K, called s, is one 
of the sides of Ti and call j3 the angle of T 2 opposite to di . After a rigid movement we can 
assume that the vertex Vi corresponding to /3 is placed at the origin and that K is contained 
in the upper half-plane. We can also assume that V 2 is placed at the point (a, 0) with a > 0, 
being the side V 1 V 2 opposite to the shortest side of K. Define now b = | Vi V 4 I sin(/3). Then, 
we have that V 4 is placed at (cot(/3)6, 6 ). Let us notice that /? is away from vr, since /3 < v/’m- 
Moreover, it is also away from 0 as one can see by means of the law of sines and taking into 








THE MINIMAL ANGLE CONDITION FOR QUADRILATERAL FINITE ELEMENTS 


10 


account that di (the side of T 2 opposite to /3) is comparable to the largest side of T 2 (due to the 
fact that 1 ^ 2 ! < -^|di| and recalling that the diameter of K agrees with the length of the longest 

diagonal). Then the linear mapping L associated to the matrix B = (^ cot(/3)\ pgj.£Qpj^g 


desired transformation L{K{a, b,d,b)) = K while it fulfills the requirements of Definition [2T] as 
||i?||, ||i?“^|| < C (with C depending on ipM,N). In particular k{B) < . 2,p,\ ■ On the other 

hand, calling L(a, b) = V 3 we observe that (A2) holds since L preserves lengths (up to constants 
depending on ||i?||, || and V 3 V 4 = L{s). On the other hand, since Ti verihes MAC{'iPm) 

and s is the shortest side of Ti then the angle of Ti placed at the common vertex of di and s is 
away from 0 and vr. Therefore {D2) holds thanks to Lemma l3.ll The theorem follows. □ 

From now on (see Lemma l2.2l) we assume that any element verifying the RDP is of the kind 
K{a,b,d,b) under [A2,D2] (equiv. [A1,D2]). In [3] it is proved that the RDP is sufficient to 
get optimal order error estimates in IT^’^ for Q\ whenever 1 < p < 3. In the last section we 
give a counterexample showing in particular that this result does not hold for k >2. 

The next result, borrowed from [3], help us to shorten our exposition playing also a role in 
the construction of a counterexample. 


Lemma 3.3. Let K = K{a,b,d,b) a convex quadrilateral. Assume [/S.2,D2\ (equiv. [Al,D2]), 
then for any 1 < p < 3 


max 


6P-1’ aP-i 


Ip<C- 


ip-i 


with C a constant depending only on those given in [A2, D2] and p. 


(3.14) 


Proof. See the proof of Lemma 3.5 in [3l pag. 140] together with eqs. (15) and (16) in that 
paper. In the mentioned lemma it is precisely the expression (|3.14p what is proved. K{a, b, d, b) 
and Ip have the same meaning that in the present work (see also [3l pag. 136], where the invoked 
hypotheses (H1),(H2),(H3) and (H4) are introduced and derived from the RDP condition.) □ 


Remark 3.1. Although it is not written here we know from [3] that the constant C in (I3.14[) 
may behave like 1/(3 —p). In order to get an uniform bound for 3 < p it is necessary to restrict 
the class of the underlying reference elements K = K{a,b,d,b). Later we show that (I3.14p holds 
for any 1 < p if we work with the family K{a, b, d, b) associated to the DAC. 

3.2. The regularity condition h/p < a and the family K{a, b, d, b). For dealing with regular 
elements we need to introduce a new geometrical condition associated to the class K{a, b, d, b), 

Condition (D3): a ^ b. (3.15) 

c 

Theorem 3.2. Let K be a general convex quadrilateral. Then K is regular (in the sense of 
(|1.4p ) if and only it is equivalent to some K{a,b,d,b) under [A2, D2,D3] (equiv. [A1,D2, D3]/. 

Proof. Thanks to Remark o we know that elements satisfying the regularity condition (|1.4p 
satisfy also the RDP. Therefore from Theorem 13.II we see that K can be mapped, by means of 
an affine transformation LX = BX + P (see Definition 12.ip into an element K(a,b,d,b). This 
reference element should be regular since ||R||, ||R~^|| < C and now it is easy to see that (Al) 
implies {D3). To prove the other implication take an element K{a,b,d,b) under [A2,D2,D3] 
(equiv. [Al, D2, D3]). Thanks to {D3) we have that T{a, b) is a regular triangle. Now it is easy 
to see that this together with [Al, D2] implies that K{a, b, d, b) verifies (11.41) . As a consequence, 
any element of the form L{K{a, b, d, b)) is regular for any affine mapping of the kind considered 
in Definition 12.11 □ 
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3.3. The DAC and the family K{a, b, d, b). As it is mentioned before DAC implies the RDP, 
and as consequence Theorem 13.11 says that any element under the DAC can be mapped into 
an element K{a,b,d,b) for which [A1,Z12] holds. Nevertheless, the following result, partially 

borrowed from [3], states that actually we may assume ^ < 1, and this, as we show later, not 
only simplifies the treatment of the error but also allows to deal with the case 1 < p. 


Theorem 3.3. Let K be a general convex quadrilateral. Then K satisfies the DAC{firn,'fiM) if 
and only if it is equivalent to an element K{a,b,d,b) under [D1,D2]. 


Proof. Notice that it is always possible to select two neighboring sides li,l 2 of K such that 
the parallelogram dehned by these sides contains the element K. Call lA the common vertex of 
li,l 2 and fi the angle at Vi- After a rigid movement we may assume that Vi = (0,0) and that I 2 
lies along the x axis (with nonnegative coordinates (a, 0)). Moreover we can also assume that li 
belongs to the upper half plane. Notice that l\ is the side joining V 1 V 4 and following the proof 
of Theorem [Q take b = |/i| sin(/3) in such a way that V 4 can be written as 14 = ( 6 cot(/3), 6 ). 
Then the linear mapping L associated to the matrix B dehned in such theorem performs the 
desired transformation. Indeed, since ||i?||, ||.B~^|| < with fi away from 0 and tt (due to 

the fact that K is under the DAC) we know that L is of the class considered in Dehnition 12.11 
On the other hand, calling T(a, 6 ) = V 3 we observe that {D\) holds while thanks to the fact 
that the angle at V 3 is away from 0 and vr the same holds for the angle at (a, b) meaning that at 
least one of the remaining angles of the triangle of vertices ( 0 , 6 ), (a, 6 ), (a, 0 ) does not approach 
zero nor tt . Performing a rigid movement if necessary we may assume that this is the one at 
(0,6) and hence {D2) follows. Reciprocally, assume that K{a,b,d,b) verihes [D1,D2] and it is 
equivalent to K. Notice that the maximal and minimal angle of K{a,b,d,b) are away from 0 
and TT (in terms of the constants given by [D1,D2]). Indeed, since at Vi we have a right angle 
we only need to check the remaining vertices. The angle at I /4 is bounded above by 7 r /2 due to 
(Dl) and below by a. Let us focus now on the angle at vertex V 3 . It should be bounded below 
by 7t/2 due to {Dl). On the other hand, it can not approach 0 due to {D2). Finally, the angle 
at V 2 is greater than a and also bounded above by 7 r/ 2 . The proof concludes by using that K 
is equivalent to K{a,b,d,b) (in the sense of Dehnition 12.ip and Lemma l3.ll □ 

There is a property that can be derived from [Dl, D2]. 


Lemma 3.4. Let K{a,b,d,b) be a general element under [D1,D2], then a < Cb. 

Proof. The proof is elementary since tan a < -. □ 

An advantage of the DAC is that it simplihes the treatment of Ip. Indeed from {Dl) we get 

1 ^ 1 

(1 + x{l - 1) + y(f - l))P-i “ (I + f - ’ 

since 0 < < 1. On the other hand, calling y{x) = —{b/a){x — a) to the equation of the 

straight line joining V 2 and V 4 we have 

b - y{d) ^ d - y~^{b) ^b d_ 
b a b a ' 

and since |/| sin(Q;) < 6 — y{d) and |/| sin(Q;) < d — y~^{b) we get 


Ip<mm{aP S 6^ 


1 


\l\P 1 sin(a)P (3.16) 

As a consequence we obtain for DAC the following equivalent of Lemma 13.31 that holds for any 


1 <p. 
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Lemma 3.5. If K = K{a, b, d, b) is a general convex quadrilateral under [Dl, D2] then, for any 

1 <P, 


max 


Ir,<C, 


h 


bP-^ ’ aP-^ j ^ ~ 

with a constant C depending only on those of [D1,D2]. 


ip-i 


(3.17) 


3.4. The mac and the family K{a,b,d,b). We finish this section with a characterization of 
the elements under the minimal angle condition. A direct consequence of previous subsections 
an Remark O is the following 


Theorem 3.4. Let K be a general convex quadrilateral. Then K satisfies the mac if and only 
K is equivalent to some K{a,b,d,b) for which holds either [D1,D2] or [Al, 112,113] (equiv. 
[A2,D2,I13];. 

Proof. Follows straightforwardly from Theorems 13.21 and 13.31 together with Remark 11.41 □ 


4. Triangles under the MAC inside K{a,b,d,b). 

Let K = K{a,b,d,b) be a general convex quadrilateral and consider the associated triangle 
T = T{a,b). Let us recall that Mij (resp. Mf^) are the interpolation nodes of Qk on K (resp. 
Ilfc on T). We are interested, loosely speaking, in the problem of finding for each Mij a close 
enough M^. Notice that in general Mij does not agree with Mfj, except for i = 0 or j = 0. 
For any other node Mij (i.e. for i 0 ^ j) we consider a suitable triangle having Mij as one of 
its vertices and with the remaining vertices belonging to the set of (edge) interpolation nodes 
of Ilfc. We choose it in the following way: if Mij is an edge node on the top of K{a, b, a, b) (i.e. 
i = k,l < j < k) Yje consider the triangle = A(MfcjM^M^^-^) if Mij is an edge node on 
the right edge (i.e. j = k,l < i < k) we chose and finally if Mij is 

interior (i.e. 1 < i,j < A: — 1) we define a triangle Tij = A^MijM^jM^) (see Figure [2|). The 



Figure 2 . Representation of T32, r23 and r2i in a Q 3 element left: T32, center: 

T 23 , right: T 21 . 

geometry of these triangles are important in the sequel. In particular notice that Tkj and Tij^ 
are similar to the triangle A(V 2 V 3 V 4 ) therefore we have immediately 

Lemma 4.1. Let K{a, b, a, b) be a general convex quadrilateral under [Al, D2] (equiv. [A2, D2], 
then for any T = T^j (resp. T = Tik) defined above we have that the side M^.jM'^Q (resp. 

is comparable to I = V3V4 and the angle of T at M^q (resp. Mf^_J is the angle a 
of condition {D2). In particular T verifies the MAC. 
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For interior nodes we have the following 

Lemma 4.2. Let K = K{a,b,d,b) be a convex quadrilateral which satisfy either [D1,D2] or 
[/S.1,D2,D2)] (equiv. [/S.2,D2,D2>]). Ifl<i,j<k — 1 then 

(a) _ _ 

\MiQMij\ ~ a and |MojMjj| ~ b 
c c 

(in particular |MjoMojl ~ hj. 

c 

(b) Oij is bounded away from 0 and it where aij is the angle between and MiqMqj. 

In particular for 1 < i,j < k — 1, any triangle T = A{Mio, Moj,Mij) verifies the MAC. 


Proof. 


(a) To prove that the measure of the segment is comparable to a it is sufficient to 

prove that the measure of the segment is comparable to a since these segments 

are mutually proportional. 

For 0 < y = i/k < 1 have 

iMioMjfcp = \\FK{l,y) - FK{0,y)f = {a{l - y) + dy,{b - b)y) 

therefore 

0^(1 - yf < iMioMjfcp < a^(l - yf + 2ady{l - y) + y‘^\l\‘^. 

Using (Al) (or (T*l)) and that I is comparable to the shortest side the statement is 
proved. Similarly, to the appropriate 0 < x < 1 we have 


\MojMkjf = \\Fk{x,1) - F'k{x, 0)\\ = {{d - a)x,b{l - x) + bx) 

therefore for a suitable constant C, we get 
6^(1 — x)^ < < 2 (a — a)^x^ + 6 ^ + (6 — 6 )^x^ <C'(a^ + 5^) (4.18) 

and the proof concludes by using {D3) in one case or Lemma 13.41 in the other. Now we 

immediately have, by using Al or Dl, that |MjoMoj| ~ h. 

c 

(b) Calling pL the matrix with rows wi = Mjo — Mik and W 2 = Mjo — M^j we see that 

1 ||rci||||t(; 2 || 

sinaij |det/r| 

Thanks to the previous item we know that the numerator can be bounded in terms of 
ab. We claim that 0 < ab < C\ det /r|. Indeed, a direct calculation gives for y = |: and 

I det fi\ = aby 


1 + --1 h /+ 7-1 


and since 1 < i < A; — 1 all we need to show is that term inside the modulus m = 
1 + (I — l) y + — 1^ X stays away from zero. Now, if | — 1 > 0 then m > 1 — y and 

we are done. On the other hand, if ^ — 1 < 0 then we write 


m > 


+ - - 1 y + -(1 - y ) > -(1 - y ) 


where the first inequality follows taking x = 1 and the second one by (j2.8p . Using Lemma 
12.51 (see also Remark EH) the proof is complete. □ 
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Lemma 5.1. Let K{a,b,a,b) be a convex quadrilateral and assume (Al). LetT = T{a,b), then 
for any polynomial q G Pfc, there exists a constant depending only on k and on the C given in 
(Al) such that 


\^\\o,p,K(a,b,d,b) — ^lkllo,p,T- 


(5.1) 


Proof. The proof is standard. Let us introduce an small rectangle Kg and a large rectangle 
Ki as follows 

Kg := K{a/2, bl2^ a/2, 6 / 2 ) <ZT G K{a,b, a, 6) C Ki := K{a,b,a,b) 
where a = max{a, a},b = max{b, b}. All we need is to show that 


q\\o,p,Ki < C\\q\\o^p^Ks- 


( 5 . 2 ) 


Thanks to (Al) we have that the quotients ^ are bounded in terms of a generic constant 
C. For the sake of clarity we rename this time the constant and write C. Consider now the 
reference sets 


kc = K 


1 1 1 1 \ 


ck = K{l,l,l,l). 


Using equivalence of norms in the finite dimensional space we get 


Iloilo,p,A — ^II^IIo,p,a:c’ 

for any g G P^ and where C depends only on k and C. Now (15.2p follows by changing variables 
with a linear map L : K ^ Ki taking into account that for such an L, L{Kq) C Kg. □ 

Write K = K{a, b, d, b) and let 11^ be the Lagrange interpolation operator of order k on the 
triangle T = T(a, b) and let p > 1 then we can write 


\u Qk'^\l,p,K ^ \'^ ^k'^\l,p,K T Infclt Qk'^\l,p,K■ 

Since — QkU belongs to the Qk quadrilateral finite element space and vanishes at Mqj 
and Mjo for all 0 < < k, it follows that 

(IlfcU - Qku){X) = ^ (IlfcU - u){Mij)(j)ij{X) 
i,j¥=o 

where (fij is the basis function associated to Mij. Therefore 

I'll Q ku\i^p^K ^ |u Lt/jul ^ ^ |(IIfcrt rt) (Afjj) 11 11 p PC. (^A) 

i,j¥=o 

Taking into account that T verifies the MAC (actually MAC{'k/2)) we have [3 [TIE] that 

'R|fc+l,p,T) (^-4) 

and 

\u - IlfcuIgp^T < (5.5) 

The next lemma extends this approximation result to K. 

Lemma 5.2. Let K = K{a,b,d,b) be a convex quadrilateral and assume (Al). Let T = T{a,b) 
and UkU the P^ Lagrange interpolation operator on T. Then for any 1 < p 


\u LT/jitlpppc ^ Ch Ittlfc^i^p^pC' 


( 5 . 6 ) 
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(5.7) 


/ D^u = / D^VkU {\a\ < k). 

Jk Jk 

Since K is convex, we have 

1 ^ 'k^k'^\l,p,K — Ch llllfc+ljPjA) 

as one can see by applying repeatedly the Poincare inequality. Writing 

^ \u 'Pkull^p^K I'PkU ^ku\l^p^K^ 

we observe that the hrst term is hne. For the second one we consider an arbitrary hrst derivative 
of VkU — UkU and call it D{VkU — Ilfcu) G Pfc-i. Using Lemma ISTTl 

\\D{VkU — Ilku)\\o^p^K < C\\D{VkU — nferi)||o,p,r < C [\VkU — u\i^p^K + \u — 'V.ku\i^p^T] 
and the lemma follows from (j5.5p and dsn. □ 

Lemma 5.3. Let K = K{a,b,d,b) be a convex quadrilateral. 

(a) Assume that 1 < p < 3 and that K satisfies [A1,D2] (equiv. [A2,L>2]J then for any 
basis function f, 

h^/P 

\f\i,P,K < 

where q is the conjugate exponent of p (the constant C may behave as see Remark 


EJl). 


(b) Assume that 1 < p < 3 and that K satisfies either [Al, 772,773] (equiv. [A2,772, 773]J 
then 

h^/P 

\^\lp,k < 

where is an internal basis function. 

(c) For any 1 <p, assume that K satisfies [771,772] then for any internal basis function 

h^/P 

(d) For any 1 <p, assume that K satisfies [771,772] then for any edge basis function 

h^/P 

Proof. Part (a) follows from Lemma 12.31 and Lemma 13.31 On the other hand, by Lemma 12.41 
and Lemma 12.51 we notice that to show (b) it is sufficient to prove that 

- Hwh . ( 5 . 8 ) 

Using that a < |Z| and \b — b\ < \l\, together with (A2) and Lemma 13.31 we have 
-^|1 - b/b\Plp < C^-^7p < C-^ < 

ap-i I IIP- jjp ^p-i p - - ^-1 

where the last inequality follows from (773). Similarly, 

^{d/a)Plp < C-^. 

6 P-1 ^ > P - qP-1 

Item (c) follows similarly to item (b) using Lemma 13.51 instead of Lemma 13.31 and Lemma 13.41 
instead of (773). Finally, the last item (d) follows straightforwardly from Lemma [2.3l and Lemma 
[331 □ 

Let us now recall the following 
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Lemma 5.4. Let T be a triangle with diameter hx and e he any of its sides. For any p > 1 we 
have 



where q is the dual exponent of p. 

Proof. See for instance m- ° 

Now we are ready to get bounds for |(rt — Ilku){Mij)\. In order to do that we consider the 
triangle Tij associated with Mij defined in Section 01 

Lemma 5.5. Let K = K(a,b,d,b) be a convex quadrilateral satisfying either [D1,D2] or 
[Al, D2, D3] (equiv. [A2, D2, D3]). For any p > 1 we consider its dual exponent q. We have, 

(a) (Edge nodes) Assume either i = k and 1 < j < A: or j = k and 1 <i <k then 

where T = Tij. 

(b) (Interior nodes) Ifl<i,j<k — 1 

nftq. 

|(n-nfcn)(M,,)|<C^ 
where T = Tij. 

Proof. 

(a) We write the case i = k and 1 < j < A: since the other one follows identically. Calling e 
the side of T = T^j given by e = Mk^M^j we get (by using Holder’s inequality, Lemma 
15.41 and the fact that (u — lij-yfiFIko) = 0) 

\{u -'Ilku){Mkj)\ < j\de{u - Ilku)\ dx 

< |e|V'/||5e(n-nfcn)||o,p,, (5.9) 


\u-Yiku\ipT + hT\u-Iiku\2^p^T > 
then 

\u - Ukuli pT + hx \u - Ilku\2^p^T > 


< 


ol/o |e| 
|T|1/P 


\\de{u 


nfc«)|lo,p,T + \de{u - nku)\;^px 


The item follows now by Lemma 14.11 that implies < C . 

(b) With the same ideas, consider {u — Ilku){MiQ) = 0 call e = Mi^Mij and use now Lemma 

K2[ □ 


Lemma 5.6. Let K be a general convex quadrilateral and 1 < i,j < k, 

(1) If 1 < p < 3 and K satisfies [Al, D2, D3] (equiv. [A2,D2,D3]), then (I5.10p holds. 

(2) If 1 < p and K satisfies [D1,D2], then (|5.10l) holds. 


I(rt Hku)[Mij^\\(j)ij\\x,K 'Ll Ch (5.10) 

where (fij is the function basis associated to Mij. 

Proof. The proof is essentially a combination of Lemmas 15.31 and 15.51 together with the error 
estimation for triangles (15.51) . recalling that each Tij satishes the maximum angle condition 
(Lemmas 14.11 and 14.21) . □ 
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6. Main Theorem 

The L'P error estimate for a general convex quadrilateral was done in [3] for k = 1 and any p. 
The argument used there works exactly in the same way for an arbitrary k. 

Theorem 6.1. Let K be an arbitrary convex quadrilateral with diameter h. For any 1 < k and 
1 < p. There exists a constant C independent of K such that 

\\u - Qfc^^|lo,p,R < Ch’^^^\u\k+i,p,K- (6.1) 

Proof. See equation (41) in Theorem 6.1 of [3], as well as Lemma 6.1 in the same paper. □ 
Now we can present our main result. 

Theorem 6.2. Let K be a convex quadrilateral with diameter h and 2 <k an integer: 

(1) If K satisfies hence (|6.2I) holds for any I <p with C = 

(2) If K satisfies mac{'iprn), hence dlJ]) holds foT oTiy 1 ^ <c 3 with C — 

\\u - QkU\\o,p,K + h\u- QkU\i^p^K ^ Ch^+^\u\k+i,p,K- (6.2) 

Proof. Since the IT estimate holds for any convex quadrilateral, it is enough to prove 

1^ Qk^h^p^x — (6.3) 

Moreover, in order to prove (1) (resp. (2)) and thanks to Theorem 13.31 freso. Theorem 13.41) 
we can assume that K = K{a,b,d,b) under [D1,D2] (resp. either [A1,H2,Z13] or [HI, 112]). 
Therefore (|6.3p follows from (15.3p combined with (j5.6p and (|5.10p . □ 

To finish we present two counterexamples. In the first one we focus on the case 1 < p < 3 
showing a collection of elements with uniform RDP parameters (actually RDP{^/ 5 , 3 / 4 : 7 r)) for 
which the constant in the interpolation error blows up. The family does not obey mac 

although all the elements are under the MAC{^). In particular, the counterexample shows 
that the estimate may fail if an angle approaches cero. For the sake of simplicity we choose 
k = 2 . 

Counterexample 6.1. (Case 1 < p < 3j For 0 < s < 1/2 take K = iL(l,s,s,2s) and 


(s,2s) 



Figure 3. Representation of the quadrilateral iF(l,s,s,2s) and its nodes as a Q 2 element. 

consider the function u{x,y) = x{x — l/2)(x — 1) which does not belong to the Q 2 space. Since 
u{Mqi) = 0 = u{Miq) for 0 < I < 2 we have 

Q 2 U = u{Mii)(f)ii + u{Mi 2 )(j)l 2 + u{M22 )4>22 + u{M 2 l)(j )21 
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and therefore 

1 dQ 2 U _ {s - 3){s + 1) d4>u f s + l d4>i2 , ^ . d<p 22 {s- 2 ) d4>2i \ 

s — 1 dy 2® dy ^ \ 2^ dy ^ dy 2^ dy ) 

Then, for a suitable constant C independent of s, we have 


dcfii 

< C 

dQ2U 

.s[ 

dcfu 

+ 

d(l)22 

+ 

d4>2i 

)] 

dy 

Ci,p,K 

dy 

0,p,K \ 

dy 

0,p,K 

dy 

0 ,p,K 

dy 

0 ,P,kJ _ 


Using item (a) of Lem,ma \5.S\ and taking into account that ~ 1 and \l\ ^ s we have 


d(pu 


dy 


< 3C 


0,p,i^ 


dQ2U 


dy 


+ s 


i/p 


0,p,K 


Assume that (insD holds for this family. In that case we would have 

< \Q2U — u\i^p^K < Ch‘^\u\3^p^K 


dQ2U 


d{Q 2 U - u) 

dy 

0,p,K 

dy 


0,p,K 

where h? and \u\^^p^K ~ ~ Consequently 


dcfii 


dy 


< Cs^/P. 


0,p,K 


On the other hand, a straightforward computation shows that 

'd(pu _ ^ \ 2^x[{s - l)y{x - y) + (1 - x)(l - 2y)] 


dy 


-oFk] {x,y) = 


Therefore 


d(j)u 


dy 


j 

J\o. 


s[l + X + (s - l)y] 

2^PxP\{s - l)y{x - y) + (1 - x)(l - 2y)|P 


o,P,K ^[o,i]2 sP i[l+ X +(s - l)y]P 1 

Let R = [0,1/8] X [1/4,3/8] C [0,1]^. It is easy to check that on R we have 
{s - l)y{x - y) + (1 - x)(l - 2y) > {s - l)y{x - y) > 0 
which together with the fact 1 + x + (s — l)y < 1 + x allow us to obtain 
d(j)ii 


dxdy 


dy 


> f r rfi#. 


sP 


-1 


(6.4) 


(6.5) 

( 6 . 6 ) 


o,p,A - h [l + x]P-i 

Since the function yP(x — y)P is bounded below by a positive constant on R and the function 
xP/(l + x)P“^ is integrable over this domain follows that 


d(fu 


dy 


> C- 


1 


0,p,K 




(6.7) 


where q is the dual exponent of p. Finally, combining with |g. ?D and taking s ^ 0 we 

are lead to a contradiction and as a consequence the error estimate can not hold with a uniform 
constant C. 


Remark 6.1. Recall that for k = 1 and 1 < p < 3 the constant in the interpolation estimate 
can be bounded in terms of the constants given in the RDF condition [3]. Actually, the interior 
node (available in k = 2) plays a fundamental role in the counterexample. This could lead the 
reader to the conclusion that removing internal nodes may help to weaken the conditions under 
which the estimate dOD holds. Regretfully this is not possible. Indeed from [4] we know that the 
accuracy of serendipity elements can be seriously deteriorated even for regular elements. The 
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reason of that is the failure of the inclusion of in the interpolation space. Our proof relies 
strongly on this property (see for instance the derivation of (|5.3I) ). 


(0.1) 



Figure 4. Representation of the quadrilateral K(l, 1, s, s) and its nodes as a Q 2 element. 

Counterexample 6.2. (Case 3 < p) Consider the family K{1, l,s,s), with ^ < s < 5/8, and 
the function u{x,y) = x{x — l/4)(3: — 3/4)(x — 3/8)(x — 1). Observe that the maximum angle 
of K = s, s) approaches tt as s ^ ^ while K verifies mac( 7 r/ 4 ) for any value of s in the 

selected range. Arguing as in previous counterexample we have 

Q2U = + u{Mi2)(t)l2 + u{M2l)(f2l + u{M22)(f'22, 

hence 

d{Q2u) ,(9012 , ,9021 , ,,,, ,9022 

Observe that u{Mii),u{Mi 2 ) and u{M 2 i) are polynomial expressions in the variable s having 
1/2 as a single root. Therefore we can write 

u{Mij) = {s- l/2)qij{s) 

for each {i,j) = {( 1 ) I); ( 1 ; 2), (2,1)} where qij is a polynomial. On the other hand\u{M 22 )\ > 

C>0if^<s< 5/8. Therefore 


9022 


dy 


< C 


0,p,K 


dQ 2 U 


dy 


+ (s —1/2) ^ 


dcfij 


dy 

0,p,K 


If the error estimates holds then 


dQ2U 


dQ 2 U — u 

dy 

0,p,K 

dy 


0,p,K 


< \Q2U - u|l,p < C\u\s,p ,Kj 


since ^ ~ 1 o.'iT'd as a consequence 


dQ 2 U 


dy 


< C. 


0,p,K 


( 6 . 8 ) 


(6.9) 


On the other hand for 1/2 < s < 5/8 we readily notice that sin(Q;) ~ 1/2)- Then combining 
this with (I3.16|) and Lemma lil.dl we get 

dc/ij 


E 


dy 


< c -- -r- 

0,p,K ('® “ 1/2)^/*? 


( 6 . 10 ) 

















































THE MINIMAL ANGLE CONDITION FOR QUADRILATERAL FINITE ELEMENTS 


20 


where q is the dual exponent to p. Finally, 116.10\} combined with id. 8\} with 116.9\) give us 

d4>22 


dy 


< C 


( 6 . 11 ) 


OjPjA 


for some positive constant. However, a straightforward calculation yields 


f <9022 
V dy 


o Fk {x,y) = 


2x [(s - l)y{x -y) + (2x - l)(4y - 1)] 


hence 


0022 


dy 


1 + {s - l){x + y) 

^ {2x [(s - l)y{x -y) + {2x - l)(4y - 1)])^ 


0,p,K JO Jo 


dxdy. 


(1 + (s - l){x + y))P-^ 

Let T be the triangle with vertices (3/4,3/4), (3/4,1) and (1,1). It is easy to check that 

0022 „ f 1 


dy 


> C 


0,p,K 


I 


T (1 + (s - ^){x + y))P 


— dxdy, 


and integrating explicitly for p > 3 we get 


0022 


and hence 


dy 

0022 


> C 


{2s - lf-P/2 + (3s - 1)3 -p/24-p - {7s - 3f-P/4^-P 


dy 


o,p,K {s-iy{2-p){3-p) 

—>■00 if s ^ 1/2. Since this fact contradicts i6.11\) we conclude that the 


0,p,K 


error estimate does not hold. The case p = 3 follows similarly. 
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